London penetration depth and superfiuid density in single crystals of Fe(Te,Se) and 

Fe(Te,S) superconductors 
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The in-plane London penetration depth, A(T), was measured in single crystals of the iron- 
chalcogenide superconductors Fei.o3(Teo.63Seo.37) and Fei.o6(Teo.8sSo.i4) by using a radio-frequency 
tunnel diode resonator. As is also the case for the iron-pnictides, these iron-chalcogenides exhibit 
a nearly quadratic temperature variation of A(T) at low temperatures. The absolute value of the 
penetration depth in the T — > limit was determined for Fei.o3(Teo.63Seo.37) by using an Al coating 
technique, giving A(0) w 560 ± 20 nm. The superfiuid density p s (T) = A 2 (0)/A 2 (T) was fitted with 
a self-consistent two-gap 7— model. While two different gaps are needed to describe the full-range 
temperature variation of ps{T), a non-exponential behavior at low temperatures requires additional 
factors, such as scattering and/or significant gap anisotropy. 



The majority of recently found iron-based supercon- 
ductors are pnictides [l[. The only exception so far 
is binary iron-chalcogenide FeSe [|J that becomes su- 
perconducting with the excess Fe occupying interstitial 
sites of the (Te,Se) layers In these materials, gen- 
erally referred to as "11" compounds, Fe forms square 
planar sheets whereas Se ions form distorted tctrahedra 
surrounding the Fe ions, which is similar to the struc- 
ture of the Fe-As pnictides. The electronic structure is 
also similar to pnictides. For "11" system it has been 
suggested both theoretically Q and experimentally [H[ 
that superconductivity could be magnetically mediated. 
Furthermore, the series of iron-chalcogenides from FeS 
through FeTe was theoretically explored within the spin- 
fluctuation picture, concluding that doped FeTe could 
exhibit the strongest superconductivity [J] . The systems 
over which the doping is most controlled are FeTei_ x Se x 
@] and FeTci_ x S x 0. So far the highest T c w 15 K is 
reported for the Fe(Te,Se) system The connection 

between superconductivity and magnetism in the "11" 
system has been demonstrated by the observation of the 
antiferromagnetic order in Fei+j,Te Q and a spin reso- 
nance in Fei+ y (Teo.6Seo.4) @. Fe(Te,S) is a supercon- 
ductor with T c rj 8.8 K and its comprehensive character- 
ization is described in Ref.Hol. 

The "11" system exhibits many interesting phenom- 
ena. The transition temperature can be enhanced up to 
37 K by applying modest pressures ll|, which is com- 
parable to the T c of iron-arsenide superconductors. The 
connection between T c and the pressure has been sug- 
gested to come from the enhancement of spin fluctuations 
12| and from the modulation of electronic properties due 



ence peak in NMR measurements on polycrystalline FeSe 
suggests unconventional superconductivity [13[ , while the 
power-law temperature dependence of the spin-relaxation 
rate, 1/T\ ~ T 3 , could be reconciled with both a nodal 
gap or a fully-gapped s± state. Muon spin rotation study 
of the penetration depth in FeSe^ was consistent with ei- 
ther anisotropic s-wave or a two-gap extended s-wave 
pairing 14[. Thermal conductivity measurements con- 
cluded multigap nodeless superconductivity in polycrys- 
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to evolution of the inter-layer Se-Fe-Se separations 
Several experimental works explore pairing mechanism 
of "11" compounds [l^ - llsT ]. The absence of a coher- 



tallinc FeSe! 

In this work, we present an experimental study of the 
London penetration depth, A(T), in single crystals of 

Fei.o3(Teo.63Se .37) and Fei.oeCIfeo.ssSo.u). We found 
that at low temperatures AA(T) cx T n with n ~ 2.1 
for Fe(Te,Se) and n ks 1.8 for Fe(Te,S). The absolute 
value of A(0) ~ 560 nm was determined for Fe(Te,Se) 
by measuring the total X(T) of the sample coated with 
a thin Al film [lij]. The in-plane superfiuid density 
p s (T) = A 2 (0)/A 2 (T) was analyzed in the framework of 
a self-consistent two-gap 7— model [l7j |. 

Single crystals of Fe(Te,Se) were synthesized using a 
flux method. Mixed powders of the Fe(Teo.6Seo.4) com- 
positions were sealed in evacuated quartz tubes. The 
sealed ampoule was slowly heated up to 930 °C and 
slowly cooled down to 400 °C at a rate of 3 °C/hr before 
the furnace was shut down. Single crystals with centime- 
ter dimensions can easily be obtained with this method 
and are shown to be the pure a-phase with the PA/nmm 
space group by x-ray diffraction Q. The actual concen- 
trations were analyzed using an energy dispersive x-ray 
spectrometer (EDXS). The measured composition for the 
samples used in this work is Fei.o3(Teo.63Seo.37). Single 
crystals of Fe(Te,S) were grown from TeS self flux using 
a high temperature flux method. Elemental Fe, Te and 
S were sealed in quartz tubes under a partial argon at- 
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FIG. 1: (Color online) Main panel: Variation of the London 
penetration depth, AA(T), for three Fe(Te,Se) samples and 
one Fe(Te,S) sample. Inset, AA(T) for Fe(Te,Se) #1 in the 
low temperature range up to T c /3 shown along with the fitting 
curves assuming exponential or power-low behavior. 
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FIG. 2: (Color online) AA plotted vs. (T/T c ) 2 for three 
Fe(Te,Se) crystals in the main panel and Fe(Te,S) crystal in 
the inset in the temperature range up to T c /3. The curves 
for Fe(Te,Se)#l and #2 are shifted vertically for clarity by 
50 nm and 25 nm, respectively. 



mosphcre. The scaled ampoule was heated to a soaking 
temperature of 430-450 °C for 24h, followed by a rapid 
heating to the growth temperature at 850 °C and then 
slowly cooled to 820 °C. The excess flux was removed 
from the crystals by decanting. The end composition of 
the sample used to measure the penetration depth was 
Fei. O6 (Teo. 88 S o .i4)0. 

The in-plane London penetration depth, A(T), was 
measured by using a self-oscillating tunnel-diode res- 
onator (TDR) [i3,[li^23. The sample was mounted on a 
sapphire rod which was inserted into a tank-circuit induc- 
tor. The weak ac magnetic field H ac ~ 20 mOc produced 
by the coil is much smaller than the lower critical field 
H c i < 100 Oe, so the sample was in the Meissner state 
and its magnetic response was determined by the London 
penetration depth. To probe the ab— plane supercurrent 
response, the sample was placed with its crystallographic 
c— axis along H ac . The shift of the resonance frequency, 
A/ = f(T) — fo, is measured to obtain the total magnetic 
susceptibility X (T) via Af = -G4n X (T). Here / ~ 14 
MHz is the resonance frequency of an empty resonator, 
G = foVs/2V c (l — N) is the calibration factor that de- 
pends on the demagnetization factor N, sample volume 
V s and coil volume V c . The calibration factor is deter- 
mined for each sample by measuring the full frequency 
change resulting from physically pulling the sample out 
of the coil at the lowest temperature. In the Meissner 
state the magnetic susceptibility, 47rx, can be written 
in terms of A and the effective sample dimension R as: 
47rx = (X/R) tanh(i?/A) — 1, from which A can be ob- 
tained [19j . 

The main panel in Fig.Q] shows the full-temperature 
range penetration depth for Fe(Te,Se) and Fe(Te,S) su- 
perconductors. While the Fe(Tc,Se) samples show a 



fairly sharp superconducting transition with AT C < 1.5 
K reflecting the high quality of the single crystals, the 
Fe(Te,S) crystal shows a broader transition. The "maxi- 
mum slope", T^ lope , determined by taking the maximum 
of the derivative dAX(T)/dT gives T c s1o p c « 12.0 K for 
Fe(Te,Se) and T c s1 °p° w 7.5 K for Fe(Te,S). As for the 
onset values, T° nsct 13 K for Fe(Te,Se) and T c onsot w 8 
K for Fe(Te,S). The low-temperature variation of A(T) 
is examined in the inset in Fig. Q] The dashed line rep- 
resents the best fit to a standard s-wave BCS function, 
AA(T) = A(0)^7rA /2Texp(-A /T), with A(0) and A 
being free fitting parameters. The experimental data do 
not show any indication of saturation down to 0.04T C 
and the fit is not adequate. Also obtained from the fit is 
Ao = 0.5T C , which is impossible in a single-gap scenario, 
hence ruling out conventional s-wavc BCS superconduc- 
tivity. We will come back to a multi-gap s-wavc fitting 
later in the discussion of the supcrfluid density. On the 
other hand, fitting with the power-law, AA(T) oc AT n , 
n = 2.1 ± 0.01, produces excellent agreement with the 
data. 

In order to examine how close the overall power-law 
variation is to quadratic, we plot AA versus (T/T c ) 2 for 
Fe(Te,Se) in the main panel of Fig. [2] and for Fe(Te,S) in 
the inset. All samples follow the AA(T) cx T 2 behavior 
rather well. To probe how robust the power n is, we 
performed a data fit over a floating temperature range, 
from T=0 to T up , using a functional form of AA(T) = 
ciq+AT 71 . The difference between the ao term determined 
from an extrapolation from the T 2 plot in Fig. [2] and the 
power-law fit turned out to be negligible, 1.5±0.5 nm, 
and had no significant effect on the fit. The dependence 
of the other fitting parameters, n and A, on T up (selected 
in the range from T c /6 to T c /3) is summarized in Fig. [3] 
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FIG. 3: (Color online) Exponent n and pre-factor A obtained 
by fitting to AX(T) oc AT' 1 for various upper temperature 
limits shown on the x-axis. The exponents in the upper panel 
were obtained with n and A both being free parameters. In 
the lower panel, A was acquired with a fixed n. 



The upper panel of Fig. [3] shows the exponent n, which (i) 
does not depend much on the selection of the upper limit 
of the fitting range, (2) is somewhat smaller in Fe(Tc,S), 
n w 1.8, than in Fe(Te,Se) n « 2. The pre-factor A 
obtained from the fit does not depend much on the fitting 
range either. 

To calculate the supcrfluid density, we need to know 
the absolute value of the penetration depth, A(0). We 
used the technique described in Ref. A thin alu- 
minum layer was deposited using magnetron sputtering 
conducted in an argon atmosphere. The Al layer thick- 
ness, t = 100±10 nm, was determined by using an Inficon 
XTC 2 with a 6 MHz gold quartz crystal and later di- 
rectly measured by using scanning electron microscopy 
on the edge of a broken sample. By measuring the fre- 
quency shift from T <C T A1 to T > T^ 1 and converting 
it into the effective penetration depth of the coated sam- 
ple, A e ff , one can extract the full penetration depth of the 
material under study from 



A e ff = Aa 



A + AAitanh(t/AAi) 
Aai + Atanh (£/Aai) 



(1) 



where A is the unknown penetration depth to be deter- 
mined. Figure 2] shows the measured A e ff(T) that is com- 
pared to the data without Al coating. The negative offset 
of 0.05 (im accounts for the thickness of the Al layer and 
\ai(T < T C A1 ). According to Eq. (P, data plotted this 
way give the actual A(T) and its extrapolation to T = 
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FIG. 4: (Color online) Effective penetration depth in single 
crystal Fe(Te,Se) before (blue circles) and after (red triangles) 
coating with an Al layer. The curve is shifted up according 



to Eq. (TfJ and the data are extrapolated to T ■ 
fit resulting in A(0) « 560 ± 20 nm. 



using a T 



gives an estimate of A(0) ~ 560 ± 20 nm for the penetra- 
tion depth of Fe(Te,Se). More details on the method can 
be found in Ref.lla. 

The superfluid density, p s {T) = A 2 (0)/A 2 (T), shown 
in Fig.[5j exhibits a noticeable positive curvature at el- 
evated temperatures, similar to MgB 2 [ll[. This sug- 
gests a multi-gap superconductivity, which we analyze 



in the framework of the self-consistent 7— model [17]. 
Although this s— wave model cannot explain the power- 
law behavior at the lowest temperatures, it should pro- 
vide a reasonable description at elevated temperatures, 
~ T c /3 <T <T C , for which effects of gap anisotropy are 
smeared. Fitting in the temperature range from 0.45T C 
to T c , shown by a solid (red) line in Fig.[5j produces a 
good agreement with the data. To limit the number of 
the fitting parameters, the partial densities of states were 
chosen to be equal in the two bands, n\ = 0.5, and the 
first intra-band coupling parameter, Ai = 0.5, was chosen 
to produce a correct T c w 12 K assuming a Debye tem- 
perature of 400 K. The variation of Ai does not affect 
the fitting quality or relative ratios of the fitting parame- 
ters. The parameters obtained in the fit are: A2 = 0.347, 
A12 = 0.096 and 7 = 0. This result means that p s (T) at 
temperatures of the order of T c is fully described by only 
one component, determined by the band with a smaller 
gap. The existence of the larger gap and small intcrband 
coupling, Ai2, are needed, however, to maintain a high 
T c . The fit over the entire temperature range reveals a 
clear deviation from this clean exponential model at low 
temperatures. The new fitting parameters of A2 = 0.281, 
A12 = 0.117 and 7 = 0.157 are close to the previous set, 
albeit with small, but finite 7 indicating 16 % contribu- 
tion of the larger gap to the total superfluid density. The 
temperature dependent gaps obtained self-consistcntly in 
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FIG. 5: (Color online) Superfluid density p s (T/T c ) for 
Fe(Te,Se)#2 calculated with experimental AA(T) and A(0) = 
560 nm. The solid (red) line is a fit the two-gap 7— model 
from 0.45T C to T c . The dashed (blue) line is the fit over the 
full temperature range. Inset: temperature dependent super- 
conducting gaps calculated self-consistently during the fitting. 



the fitting are shown in the inset to Fig.[5J While the 
fitted positive curvature and reasonable coupling param- 
eters indicate a multi-gap nature of superconductivity 
in "11" iron-chalcogenide superconductors, the failure at 
low temperatures and apparently non-exponential behav- 
ior requires extension to the anisotropic gap and inclusion 
effects of (possibly strong) pairbrcaking [22j. 

In conclusion, a robust power-law behavior of 
the low-temperature A(T) is found in single crys- 
tals of Fei.o3(Te .63Se .37) and Fei. 6(Te .88So.i4) with 
AA(T) oc T n with n rj 2.1 and w 1.8, respectively. For 
Fe(Te,Se), the absolute value, A(0) « 560 ± 20 nm, was 
determined by the coating technique. The analysis of the 
superfluid density shows a clear signature of multi-gap 
superconductivity and a failure of the clean limit s-wave 
(including s±) pairing. 
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